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We investigate magnetic properties induced by a spinless impurity in strongly correlated electron 
systems, i.e. the Hubbard model in the spatial dimension D = 1,2, and 3. For the ID system 
exploiting the Bethe ansatz exact solution we find that the spin susceptibility and the local density 
of states in the vicinity of a spinless impurity show divergent behaviors. The results imply that the 
' induced local moment is not completely quenched at any finite temperatures. On the other hand, the 

spin lattice relaxation rate obtained by bosonization and boundary conformal field theory satisfies a 
relation analogous to the Korringa law, 1/T\T ~ x 2 . In the 2D and 3D systems, the analysis based 
| upon the antiferromagnetically correlated Fermi liquid theory reveals that the antiferromagnetic 

spin fluctuation developed in the bulk is much suppressed in the vicinity of a spinless impurity, and 
thus magnetic properties are governed by the induced local moment, which leads to the Korringa 
law of l/Ti. 
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I. INTRODUCTION 



i 

Recently magnetic properties induced by spinless impurities in correlated electron systems have attracted much 
interest JM Especially, to probe antiferromagnetic correlations of High-T c cuprites the substitution of Cu sites with 
non-magnetic impurities such as Zn, Al, and Li, has been studied experimentallyJjcl According to NMR experiments, it 
was found that the substitution with spinless impurities induces local moments in the vicinity of impurities, which also 
show Kondo-like behaviors. For instance, the spin susceptibility in the vicinity of an impurity shows the temperature 
dependence like ~ 1/(T + Tk), which implies the existence of the characteristic energy scale Tk analogous to Konda 
temperatureE Moreover the spin lattice relaxation rate l/Ti shows Korringa-like behaviors, \ jT\T cx K 2 , for T < Tk-Q 
i Here K is the Knight shift. It is noted that in the vicinity of a spinless impurity the antiferromagnetic spin correlation 

which is developed in the bulk is much suppressed, and the magnetic correlation is dominated by the induced local 
moment. From theoretical points of view, it is non-trivial how this induced local moment governs the magnetic 
properties around an impurity, suppressing the antiferromagnetic correlation. In this paper, we shall deal with this 
issue. Although the experiments are carried out for High-T c cuprates which are essentially quasi-two-dimensional 
£^ ' systems, it is expected that such effects may depend on the lattice structure and the dimensionality. Thus, we 
, consider the Hubbard models with a spinless impurity in the spatial dimension D = 1,2, and 3 to investigate how the 
C^) • dimensionality affects the induced magnetic properties. For D = 1, the effects of a spinless impurity is incorporated 
into an open boundary condition as will be explained in the next section. Thus we consider the ID Hubbard model 
with boundaries which is exactly solvable in terms of the Bethe ansatz method. We analyze the magnetic properties 
of this model using the exact solution and boundary conformal field theory. For D — 2 and 3, we derive the Korringa 
relation satisfied in the vicinity of a spinless impurity which is observed in NMR experiments. Our argument for 
D = 2, 3 is based upon the Fermi liquid theory in the presence of antiferromagnetic spin fluctuations. 

The organization of this paper is as follows. In Sec. II, the ID Hubbard model with a spinless impurity is considered. 
The spin susceptibility and the local density of states in the vicinity of an impurity are obtained based upon the Bethe 
ansatz exact solution. It is found that the induced moment is not screened completely at any finite temperatures. We 
also derive the spin-lattice relaxation rate l/Ti which satisfies a relation analogous to the Korringa law. In Sec. Ill, we 
discuss about the 2D and 3D systems exploiting the antiferromagnetically correlated Fermi liquid theory. Summary 
is given in Sec. IV. 



II. A SPINLESS IMPURITY IN THE ONE-DIMENSIONAL HUBBARD MODEL 



A. Mapping to the Hubbard model with boundaries and the Bethe ansatz exact solution 

The effects of a single impurity in one-dimensional correlated systems have been extensively studied so far.Bi If 
the interaction between fermions is repulsive, a potential scattering in the ID Hubbard model is renormalized to an 
infinite strength, eventually, cutting the system into two half-infinite chains in the low-energy scaling limit. Thus at 
sufficiently low temperatures the system can be treated as the Hubbard chain with open boundaries, of which the 
hamiltonian is given by, 
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where the last term is a boundary potential. As we will see below, the low-energy spin dynamics around the impurity 
with which we are concerned are mainly described by this model, and the interaction or hopping between the two 
half-infinite chains is a subleading irrelevant interaction which can be incorporated by perturbative calculation. n .— . 

The Bcthc ansatz exact solutions of ID correlated systems with boundaries have been studied by many authors.ETKll 
In connection with the spin dynamics in the vicinity of the boundary, an intriguing result was obtained for the 
supersymmetric t-J model by Essler.tj He obtained the divergent behavior of the boundary spin susceptibility as a 
function of a magnetic field H, i.e. Xboundary ~ l/H(lnH) 2 . It was-first predicted by de Sa and Tsvelik that such a 
Curie-like behavior is universal for integral models with boundaries .E-a Later, the similar behavior was also found for 
the Hubbard model at half-filling by Asakawa and Suzuki.E3 In the next subsection, we shall show that this divergent 
behavior holds also for the case away from half-filling with finite U. 

Here we summarize the basic equations which are relevant to the following arguments. The Bethe ansatz equations 
of the ID Hubbard model with boundaries obtained by Schulz many years ago arejill 
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U /A, and 0o,l is a potential at boundaries. TV is the total number of electrons. M is 



the total number of down spins, kj and A Q are rapidities for charge and spin degrees of freedom, respectively. In 
the following, we consider only the case of repulsive boundary potentials. Thus the above equations have real roots. 
Putting k-j = —kj, A-c = — Aq,, and taking a continuum limit, we have the integral equations for the distribution 
functions of rapidities, 
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Then the magnetization is expressed as, 
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The total energy is expressed in terms of the dressed energies, 
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where the dressed energies e c (k) and £ S (A) are determined by the integral equations, 
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If one fixes the magnetic field H, B is determined by the equilibrium condition dE/dB — 0, which is equivalent to the 
condition, s s (B) — 0. In the subsequent sections, we calculate the spin susceptibility and the local density of states 
using the above equations. 



B. Spin susceptibility 

In order to derive the spin susceptibility, we solve eqs. (@) and H for er(A) using the Wiener-Hopf method, and 
obtain the magnetization, eq.(^). 

Applying the Fourier transformation and shifting the argument, A —> A + B, we rewrite eq.(^|) into, 

f°° fl\' C°° dX' 
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The last term of the right-hand side of eq.(|l2|) is 0(1/ B 2 ) forpSpiall magnetic fields. Thus we neglect it. Then eq.fllz 
can be solved by using the standard Wiener-Hopf methodE3 The solution is expressed in terms of the following 
functions, 

G+(uj) = V2^ y ' -e^, (15) 
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/oH = / dXf (X + B)e iuX , (18) 



where Q + (u>) (Q~(u>)) is the analytic part of G~(u})fo((j) defined in the upper (lower) half plane. Fourier transforming 
eq.(|I^) and introducing the function <j + (uj) = J °° dXe tujX (j(X + B), we obtain the solution as, ct + (uj) = G + (u)Q + (uj) . 

Now we derive Q + (cu) as follows. For small magnetic fields, i.e. large B, and A > the second term of /o(A + B) 
is approximated as, 
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This driving term is essentially the same as the bulk contribution, with which we are not concerned. The first term of 
fo(X + B) gives rise an interesting boundary effect. Using e~"' a, '/2 cosh itu = ^^ i (_i)«- 1 e - 2llu l w l anc [ the Laplace 
transformation, 
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Finally, using eq. (||) , we obtain the magnetization, 
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for large B, i.e. small magnetic fields. B is related to H from the condition e s (B) = 0. From eqs.([To|) and (0), we 
have, H = Ce~* B / 2u for H < u. Here C is an constant. Then the spin susceptibility \ — dS z /dH behaves like, 

X ~ I g(lnij-) 2 + bulk temS ' (24) 

This ii-dependence is the same as that found for the half-filling case.0 The above result implies that in ID systems 
the magnetic moment induced by a non-magnetic impurity is not screened completely even at zero temperature^. This 
behavior is analogous to the underscreening multi-channel Kondo effect, as pointed out by de Sa and Tsvelik.li-J The 
leading if-dependence of eq. (^) is not altered, even if one includes irrelevant interactions such as the hopping between 
the two half-infinite chains. 

In this section, we restrict our discussion to the zero temperature case. It is expected that at finite temperatures 
the boundary spin susceptibility behaves like Xboundary ~ l/T(lnT) 2 . In order to confirm this prediction, we need to 
explore thermodynamic Bethe ansatz method in the presence of boundaries. However, in the presence of boundaries, 
the entropy can not be expressed in terms of rapidity distribution functions in the continuum limit, because of the 
presence of spurious states for vanishing rapidities, and thus the usual technique of thermodynamic Bethe ansatz 
method is not applicable. If we limit the argument to sufficiently low temperature regions, undesirable contributions 
from the spurious state around the bottom of the energy spectrum may be small, and not give rise serious errors. 
Even if we admit this approximation, it is still a cumbersome task to solve thermodynamic Bethe ansatz equations 
numerically for low temperatures. Thus here we just give a field-theoretical argument to justify the above speculation. 
According to the boundary conformal field theory, the abe*je-idivergent behavior of the spin susceptibility is due to 
the presence of a boundary entropy, Abound = T ln( y/AnR) ffl Here R is the radius of the boson field of the Gaussian 
model which is the low-energy effective theory. If the leading irrelevant interaction is the marginal operator in the 
spin degrees of freedom, J l ■ Jr, we have R ~ Rq — g/\sxT for small TEj Then, the boundary spin susceptibility 
should behave like, Xboundary ~ l/T(lnT) 2 . Thus we expect that this temperature dependence which signifies the 
presence of an unquenched local moment may realize in this system. 



C. Local density of states 



In models solvable in terms of the Bethe ansatz method, the local density of states is defined as the derivative of 
the quantum number, which parameterizes rapidities, with respect to the pseudo-energy, i.e. dlj /de(kj)SB For the 
ID Hubbard model, we can consider the density of states of holon and spinon, respectively. An interesting singular 
behavior due to the boundary appears in the spin degrees of freedom. 

The local density of states of spinon as a function of energy is given by, 
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In the absence of magnetic fields, B — > oo, the solution of eq.(H) is expressed as, 
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For A ^> 1, the first term of eq. (|26|) behaves like ~ 1/A 2 , while the second term is just the order of 0(e" 7rA / u ). Thus 
the main singular contribution comes from the former which is nothing but the boundary term. In a similar manner, 
from eq. ( [il|) we obtain the asymptotic form of e s (A) for large A, i.e. e s (A) ~ Ae~ 7rX / u , where A is a constant. Then 
from eqs.(|25D and (p6|), we have, 

Pspin(£) ~ 7ohr> (27) 



for small e. Thus the local density of states also shows the singular divergent behavior because of the presence of the 
boundary. It is noted that this result is similar to that of the underscreened multi-channel Kondo effect! 
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The important message of this and the previous subsections is that in ID correlated electron systems the localized 
moment induced by a non-magnetic impurity is not quenched at any temperatures. The inclusion of irrelevant 
interactions such as hopping between semi-infinite chains does not change the result qualitatively. It should be 
stressed that this unqucnched local moment is a particular property of the ID systems where an impurity divides the 
system into two semi-infinite chains. Such a separation of the system is not possible in higher dimensional systems. 



D. Spin lattice relaxation rate 



Here we calculate the spin lattice relaxation rate 1/Ti in the vicinity of a spinless impurity, i.e. a boundary, using 
bosonization method and boujadajy conformal field theory. The same kinds of the calculations have been done for 
Heisenberg spin chains before. Ej€3 Some parts of the following results are similar to those obtained in ref.26 and 27. 
However combining them with the results from the Bethe ansatz exact solution, we shall see some new aspects. In the 
previous subsections, it was shown that the induced moment is not screened completely at any temperatures. Then 
one might expect that 1/Ti behaves like that of an isolated spin, 1/Ti ~ T\- However, as will be seen below, this 
naive expectation is incorrect. 

According to the boundary conformal field theory, correlation functions for any operators in the vicinity of bound- 
aries are obtained by the analytic continuation of. tha. antiholomorphic part to the holomorphic part, 0(z, z) — 
OL(vt + ix)On(vt — ix) ~ Oi,(vt + ix)OL(vt — ix)E5t3 Following the standarAjt^hpique, we have the asymptotic 
behaviors of the spin-spin correlation function in the presence of the boundary,! 
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Here v s and v c are the velocities of spinon and holon, respectively. K c is the Luttinger liquid parameter in the charge 
sector, and 1/2 < K c < 1. K s = 1 because of the SU(2) symmetry of the spin sector. In the vicinity of the boundary, 
i.e. x, y, \x — y\ -C v s t, the staggered part (the second term) of eq.(p8|) is less relevant in comparison with the uniform 
part (the first term) . Thus in contrast to the bulk behavior, the antiferromagnetic spin fluctuation is much suppressed, 
and the uniform part gives the dominant contribution to 1/Ti near the boundary. Fourier transforming eq.(p8[) J^jrEHl 
we obtain, up to logarithmic corrections, 

i-=lim^Im X ( M »4 + 0(T^). (29) 
q,q' 



Here we have assumed that the hyperfme coupling constant is independent of q, and omitted it. C is a temperature- 
independent constant. In the case that v s is-a constant, the above result is equivalent to that obtained by Brunei 
et al. for Luttinger liquids with boundaries.E3 The spinon velocity v s is related to the spin susceptibility obtained 
before, l/v s = Xbulk + Xboundary/T- As claimed in the previous subsections, Xboundary should show enhanced local 
correlations like Xboundary ~ l/T(lnT) 2 . Thus near the boundary, 

TjTTj; ~ (Xboundary) 2 - (30) 

Surprisingly, this relation is analogous to the Korringa relation. However it is noted that in contrast to the conventional 
Korringa law, the right-hand side of eq. (j3TJ) shows strong temperature dependence. As mentioned in the previous 
subsections, the induced local moment is not quenched completely. In spite of such an unscreened character of the 
moment, a la Korringa relation holds in the vicinity of a spinless impurity. 



III. A SPINLESS IMPURITY IN THE 2D AND 3D HUBBARD MODELS 



In this section, we discuss the local magnetic properties caused by a spinless impurity in the 2D and 3D Hubbard 
model in the presence of bulk antiferromagnetic fluctuations, i.e. very close to the half-filling. The main purpose of 
this section is to derive the Korringa gelation satisfied at the nearest neighbor of the impurity site, which is observed 
in the NMR experiment for cuprates.Q The model hamiltonian is given by, 
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Here E k — — 2t^^ =1 cosfc a (D — 2 or 3), and the last term represents a spinless impurity localized at site 0. We 
restrict our discussion to the case of the square lattice in 2D and the cubic lattice in 3D. Because of the presence of 
an impurity, correlation functions are non-local. In the case of U = 0, the single particle Green's function is given by, 
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where fi is a chemical potential. For U =/= 0, the single particle Green's function is obtained by solving the equation, 

2j[(ie n + A 1 - £k)Skk" - S fc fe''(e„) - E Q }G k "k'(£n) = &kk>- (33) 
k" 

The self-energy T,kk'(s) may be obtained by perturbative calculation in terms of U. However, in the following 
qualitative argument we do not need the explicit expression of Gkk'{sn)- 

Before discussing about magnetic properties, it is useful to sketch the spatial dependence of the density of states 
in the vicinity of an impurity. The density of states at the Fermi level is given by, 

p(x,x') = -(lA)]TlmG«,(0)e^e-* fcV (34) 

kk' 

To simplify the calculation, we consider the strong limit of an impurity potential, i.e. Eq 3> U,t. The following 
argument do not change qualitatively even in the case of a finite Eq . Using eq. ( [32] ) , we obtain the density of states at 
the Fermi level for the non-interacting system, U = 0, 
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Now 
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with J n (x), the Bessel function, and x = (ni,ri2) for D = 2 and x = (ni, n2,na) for D — 3. Note that if x or x' 
is the impurity site, the density of states vanishes, po(x, 0) = po(0, x') = 0. We can easily show that if the electron 
density is close to the half-filling, i.e. \(i\/t -C 1, then N x (p,) ~ 0((p/t) 2 ) for the site x on the sublattice which 
includes the nearest neighbor site of the impurity, x n . n . (denoted by A-sublattice) , and N x {n) ~ No(fi) for the site 



x on the sublattice which includes the impurity site (£?-sublattice). Thus from eq.(35) we immediately see that the 
local density of states around the impurity site shows strong spatial modulation similar to the Friedel oscillation. The 
period of the oscillation is ~ 1//cf which is close to the half-filling value in this case. If x and x 1 belong to A-sublattice, 
the local density of states is nearly equal to that of bulk systems, pq(x,x') ~ Pod^l — * oo, \x'\ — > oo). On the other 
hand, if x ^ and x' ^ belong to B-sublattice, po{x, x') ~ 0((p,/t) 2 ). Since the density of states at the Fermi level 
is not renormalized by electron-electron interaction, this Friedel oscillation occurs for U ^ 0. This observation leads 
to an important implication for local magnetic properties on the nearest neighbor site of the impurity, x n . n ,. Because 
of the Friedel oscillation and the bipartite lattice structure, the local density of states on all the sites surrounding the 
site x n . n . is much suppressed provided that the impurity potential is sufficiently strong. Thus the spin on the site 
x n .n. is less screened than spins of electrons in the bulk. As a result, the local spin susceptibility on the site x n . n . is 
strongly enhanced. 

Now we consider the spin lattice relaxation rate 1 jT\ in the vicinity of a spinless impurity. For simplicity, we assume 
that the hyperfine coupling constant does not depend on q. We apply the general argument from the Fermi liquid 



theory to the case with a single spinless impurityHo Then 1/Tj at the site Xi is given by, up to constant factors, 
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where A(q,k + q',k) is a three point vertex function. The diagrammatic expression of eq.(37) which is the w-linear 
term of Imx(q, q' , to) is shown in FIG.l. The detail derivation of this formula is given in ref.33 and 34. In the presence 
of strong antiferromagnetic fluctuations, it is plausible to assume that A(q, k + q', k) depends mainly on q and q' . Thus 
in the following we discard the fc-dependence of A(g, k + q' ,k). It is useful to rewrite eq.(^7|) in terms of quantities in 
the coordinate space, 
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In the case that the site Xi is far from the impurity, i.e. |a^| 3> a, where a is a lattice constant, p(xi + s, Xi—t) — > p(s+t), 
A(xi,Xi + s) — > A(— s), and thus the above expression is reduced to the usual formula of 1/Ti in bulk systems. It is 
also noted that if xi — 0, \JT\T vanishes, since A(0, s) includes G(0, Xj) which vanishes as mentioned above. Here we 
are concerned with the case that the site Xi is the nearest neighbor of the impurity site, Xi = x n , n .. 

To proceed further, we use a phenomenological expression for A{q,q'). We assume that the three point vertex 
function A(q, q') consists of the part which is strongly enhanced by antiferromagnetic spin fluctuation and the local 
part which depends on q and q' weakly, 



A(q, q') ~ A AF (q)5 q<q < + A loc (q, q'). 



(40) 



The antiferromanetically correlated part Aaf(q) bas a strong peak at q — Q, the staggered vector, and then approx- 
imated as, ReA^p('?) ~ Rex(g ~ Q) = + (^af(<Z — Q)) 2 )- Here we used the phenomenological expression 
for x(q ~ Q)£a As mentioned above, A\ oc (q,q') is enhanced by local magnetic correlations at the site x n . n ., i-e. 
Aioc(g,<7') ~ A\ oc e- iqx * " e lq ' x * •»■ Then, eq.@ is rewritten into, 
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where m — 2 for 2D systems and m = 3 for 3D systems. Since X (Q) ~ (a F ) 2 @ the first term of eq. (fid]) , which is 
the antiferromagnetically correlated part, is much suppressed compared to the second term, i.e. the local correlation 
part. Thus we obtain, 



1 



(ReA loc ) 2 



(42) 



Here we neglect all factors which are not enhanced by electron correlation. On the other hand, the local spin 
susceptibility at a; n n . is approximately given by, xioc ~ ReAi oc . Thus eq. (|42|) establishes the Korringa relation 
satisfied at the nearest-neighbor site of the impurity. As mentioned before, this relation is actually observed in NMR 
experiments.u 




k+q 2 k'+q 3 



FIG. 1. diagram of tj-linear term of lmx{qi, Q4, w). The shaded part is the three point vertex. 
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IV. SUMMARY AND DISCUSSION 



We have discussed some magnetic properties analogous to the Kondo effect induced by a spinless impurity in 
strongly correlated electrons systems. In the ID system, we have shown that the spin susceptibility and the local 
density of states near the impurity indicate divergent behaviors implying the presence of an unquenched local moment 
at any temperatures. We have also obtained the Korringa-like relation between the spin lattice relaxation rate and 
the local spin susceptibility, 1/T\T ~ (xboundary) 2 - In the 2D and 3D systems, the antiferromagnetically correlated 
Fermi liquid theory has been applied. It has been shown that magnetic properties in the vicinity of a spinless impurity 
are dominated by the induced moment rather than the antiferromagnetic spin fluctuation developed in the bulk, and 
that the Korringa relation holds at the near neighbor site of the impurity. 

The results obtained for the ID system have an interesting implication to higher dimensional systems. Suppose 
a semi-infinite 2D Hubbard model with a boundary line, which is regarded as the coupled semi-infinite Hubbard 
chains. According to the results obtained in Sec. II, it is expected that at some finite temperatures the ID-like 
strong spin correlations occurs in the vicinity of the boundary line leading to strongly enhanced density of states 
near the boundary. Such an enhanced electron correlation and one dimensionality of the boundary line may give rise 
strong fluctuations toward some surface phase transition. For instance, if there exists a pairing interaction in the bulk 
system, the paring correlation may be enhanced near the boundary leading to higher transition temperature than the 
bulk superconductivity. Actually, it is reported that Sr2Ru04 with lamellar microdomains of Ru metal shows the 
superconducting transition at the temperature higher than T c of the pure system, and that the superconductivity 
with higher-T c occurs in the vicinity of the boundary between Sr2Ru04 and Ru-metal.EJ We would like to pursue this 
possible mechanism of the enhanced transition temperature in the near future. 
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